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Abstract. We prove the convergence in any time interval of a point-particle 
approximation of the Vlasov equation by particles initially equally distributed 
for a force in with a < 1. We introduce discrete versions of the L°° 

norm and time averages of the force field. The core of the proof is to show 
that these quantities are bounded and that consequently the minimal distance 
between particles in the phase space never vanishes. 
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1 Introduction 

We are interested here by the validity of the modeling of a continuous media 
by a kinetic equation, with a density of presence in phase and speed. It others 
words, does the many particles follow the evolution given by the continuous 
media when their number is sufficiently large? This is a very general question 
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and this paper claims to give an (partial) answer only for the mean field 
approach. 

Let us be more precise. We study the evolution of N particles, centered at 
(Xi, . . . , X n ) in R d with velocities (Vi, . . . , V n ) and interacting with a central 
force F(x). The positions and velocities satisfy the following system of ODEs 

X = V u 

V = E(X t ) = ST^1 F{Xt - Xj), (LI) 
or* 

where the initial conditions (X®, V] , . . . , X®, V®) are given. The prime ex- 
ample for (jl.lj) consists in charged particles with charges a, and masses rrii, 
in which case F(x) = — x/\x\ 3 in dimension three. 

To easily derive from (jl.lj) a kinetic equation (at least formally), it is very 
convenient to assume that the particles are identical which means cti = a.j. 
Moreover we will rescale system (jl.lj) in time and space to work with quan- 
tities of order one, which means that we may assume that 

— = Vi,j. (1.2) 
rrii N 

We write now the Vlasov equation modelling the evolution of a density / of 
particles interacting with a radial force in F(x). This is a kinetic equation in 
the sense that the density depends on the position and on the velocity (and 
of course of the time) 

dtf + v- VJ + E(x) ■ V v f = 0, t eR + , x eR d , v e R d , 

F{x) = V{J p{t,y)F{x-y)dy), (13) 

P{t,y) = / f(t,x,v)dv. 



Here p is the spatial density and the initial density /° is given. 
When the number N of particles is large, it is obviously easier to study 
(or solve numerically) ()Oj) than (jl.lj) . Therefore it is a crucial point to 
determine whether (jl.Hj) can be seen as a limit of (jl.lj) . 
Remark that if (Xi, . . . , X N , Vi, . . . , V n ) is a solution of (jl.lj) . then the mea- 
sure 

1 n 

p N {t) = jjJ2 5 ( x - 5{v - 

i=i 
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is a solution of the Vlasov equation in the sense of distributions. And the 
question is whether a weak limit / of //jy solves (jl.3|) or not. If F is smooth, 
then it is indeed the case as it is proved in the book by Spohn [20]. The 
purpose of this paper is to justify this limit if 

1^)1 |Vi^)l<^ |V 2 nx)i<^, Vx^O, (1.4) 

for a < 1, which is the first rigorous proof of the limit in a case where F is 
not necessarily bounded. 

Before being more precise concerning our result, let us explain what is the 
meaning of in view of the singularity in F. Here we assume either 

that we restrict ourselves to the initial configurations for which there are 
no collisions between particles over a time interval [0, T] with a fixed T, 
independent of N. Or we assume that F is regular or regularized but that 
the norm ||F||v^i.°o may depend on N; This procedure is well presented in [T] 
and it is the usual one in numerical simulations (see [21] and [22] )• in both 
cases, we have classical solutions to (jl-H) but the only bound we may use is 

JED. 

Other possible approaches would consist in justifying that the set of initial 
configurations Xl(0), . . . , Xjy(O), Vi(0), . . . , Vjv(0) for which there is at least 
one collision, is negligible or that it is possible to define a solution (unique 
or not) to the dynamics even with collisions. 

Finally notice that the condition a < 1 is not unphysical. Indeed if F 
derives from a potential, a = 1 is the critical exponent for which repulsive 
and attractive forces seem very different. In other words, this is the point 
where the behavior of the force when two particles are very close takes all its 
importance. 



1.1 Important quantities 

The derivation of the limit requires a control on many quantities. Although 
some of them are important only at the discrete level, many were already 
used to get the existence of strong solutions to the Vlasov-Poisson equation 
(we refer to [S], [H] and (THj, [IZ] as being the closest from our method). 
The first two are quite natural and are bounds on the size of the support of 
the initial data in space and velocity, namely we introduce 

R(T) = sup Xi(t), K{T) = sup Vi{t). (1.5) 

te[0,T],i=l,...N te[0,T],i=l,...N 
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Of course R is trivially controlled by K since 



R(T)<R(0) + TK(T). (1.6) 

Now a very important and new parameter is the discrete scale of the problem 
denoted e. This quantity represents roughly the minimal distance between 
two particles or the minimal time interval which the discrete dynamics can 
see. We fix this parameter from the beginning and somehow the main part 
of our work is to show that it is indeed correct, so take 

At the initial time, we will choose our approximation so that the minimal 
distance between two particles will be of order e. 

The force term cannot be bounded at every time for the discrete dynamics 
(a quantity like F-kp^ is not bounded even in the case of free transport), but 
we can expect that its average on a short interval of time will be bounded. 
So we denote 

E(T) = sup {- [ t+£ \E(X t (s))\ds\ , (1.8) 

te[to,T-e],i=l,...,N l £ Jt J 

with for T < e 
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T 



E(T)= sup - / \E(X t (s))\ds , (1.9) 

i=l,...,N {£ Jo ) 

thus obtaining a continuous definition. Moreover we denote by E° the supre- 
mum over all % of ^(X^O))]. 

This definition comes from the following intuition. The force is big when two 
particles are close together. But if their speeds are different, they won't stay 
close a long time. So we can expect the interaction force between these two 
particles to be integrable in time even if they "collide". They just remain the 
case of two close particles with almost the same speed. To estimate the force 
created by them, we need an estimate on their number. On way of obtaining 
it is to have a bound on 

m(T) = M^W^FW' (L10) 
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The control on m requires the use of a discretized derivative of E, more 
precisely we define for any exponent (3 G ] 1, d — a [ which we note also 
satisfies (3 < 2d — 3a {(5 = 1 would be enough for short time estimates) 

AE(T) — sup sup arM»P4, (1 .U) 

with as for -E 1 , when T < e 

AE(T) = sup P 

Now, we introduce what we called the discrete infinite norm of the distribu- 
tion of the particle fiN- This quantities is the supremum over all the boxes 
of size e of the total mass they contains divided by the size of the box. That 
is, for a measure \i we denote 

11^1100,6 = 7^6 SU P {v( B oo(( x > v )> £ ))}- ( L13 ) 

(x,»)a 6 

where B^^x, v), e) is the ball of radius e centered at (x,v) for the infinite 
norm. Note that we may bound \\/in(T, -)\\oo,e by 

||/iiv(T,-)||oo, £ <(4m(T)) M . (1.14) 

All the previous quantities (except for e) will always be assumed to be 
bounded at the initial time T = uniformly in N. 

1.2 Main results 

The main point in the derivation of the Vlasov equation is to obtain a control 
on the previous quantities. We first do it for a short time as given by 

Theorem 1.1. If a < 1, there exists a time T and a constant c depending 
only on R(0), K(0), m(0) but not on N such that for some a < a' < 3 

R(T) < 2 (1 + R(0)), K(T) < 2 (1 + K(0)), m(T)<2m(0), 

E(T) < c (m(0)) 2a ' (K(O)r' (R(0)r'~ a , sup \\fi N (t, -)||oo ;£ < (8m(0)) M . 

t<T 



5 



Remark 

The constant 2 is of course only a matter of convenience. But the choice of 
a different constant is not really helpful. This result is valid only for a short 
time in essence and increasing the chosen constant for instance, increases 
only slightly the time T who in any case cannot pass a critical value. 

This theorem can be extended on any time interval 

Theorem 1.2. For any timeT > 0, there exists a function N ofR(0), K(0), 
m(0), andT and a constant C \R(0) , K (0) , m(0) , T) such that if N > N then 

R(T), K(T), m(T), E(T) < C(R(0), K(0), m(0), T). 

^From this last theorem, it is easy to deduce the main result of this paper 
which reads 

Theorem 1.3. Consider a time T and sequence /Uat(£) corresponding to so- 
lutions to (|1-H) such that R(0), K(0) and m(0) are bounded uniformly in 
N. Then any weak limit f of /i N (t) in L°°([0, T], M\R 2d )) belongs to 
L°°([0, T], L l H L°°{R 2d )), has compact support and is a solution to (J1.3j) . 

Of course the main limitation of our results is the condition a < 1 and 
the main open question is to know what happens when a > 1. However 
this condition is not only technical and new ideas will be needed to prove 
something for a > 1. 

It would also be interesting to extend our result to more complicated forces 
like the ones found in the formal derivation of [T2"j . 

The derivation of kinetic equations is an important question both for numer- 
ical and theoretical aspects. We already mentioned the works of Batt pQ, 
Spohn [20], Victory and Allen an d Wollmann [22] for Vlasov equations. 
Another interesting case concerns Boltzmann equation, for which we refer to 
the book by Cercignani, Diner and Pulvirenti [3] and the paper by Illner and 
Pulvirenti [TUj . 

On the other hand, the derivation of macroscopic equations is usually easier 
and some results are already known (although not since a very long time) 
even in cases with singularity. In particular and that is more or less the 
macroscopic equivalent of our result, the convergence of the point vortex 
method for 2 — D Euler equations was obtained by Goodman, Hou and 
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Lowengrub [7] (see also the works by Schochet [TH] and [IH])- The main 
difficulty for macroscopic systems is to control the minimal distance between 
two particles (which is not possible in the kinetic framework) as it is also 
clear in [TT] . 

Our method of proof makes full use of the characteristics and of the pro- 
cedures developed to get for the Vlasov-Poisson equation in dimension two 
and three. This method was introduced by Horst in jH] and |H] and was 
successfully used to prove the existence of strong solutions in large time in 
P^j and and at the same time by Lions and Perthame in using the 
moments (see also [5 for a slightly simpler proof and for an application 
to the asymptotic behavior of the equation). These results were extended to 
the periodic case by Batt and Rein in |3] and to the Vlasov-Poisson-Fokker- 
Planck equation by Bouchut in [2]. In particular the necessity to integrate 
in time to control the oscillations of the force also appears in the proof of 
L°° bounds for the Vlasov-Poisson-Fokker-Planck equation by Pulvirenti and 
Simeoni in jT^l • We refer to the book by Glassey jE] for a general discussion 
of the existence theory for kinetic equations. 

In the rest of the paper, C will denote a generic constant, depending maybe 
on R(0), K(0), or m(0) but not on iV or any other quantity. We first prove 
Theorem ll.il then we show a preservation of discrete L°° norms which proves 
Theorem 11.21 In the last section we explain how to deduce Theorem 11.31 



2 Proof of Theorem 11.11 

The first steps are to estimate all quantities in term of themselves. Then 
if this is done correctly it is possible to deduce bounds for them on a short 
interval of time. 



2.1 Estimate on E 

We are in fact able to estimate any time average of the force field, more 
precisely 

Lemma 2.1. For any a' with a < a' < 3, assume that 

, \ 1 
m(t Q ) < 



12 eK(t )AE(t o y 
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then the following inequality holds 

E(t ) <C{\\» N \\%*K a ' R a '- a + e d - a \\t, N \\ 00 , £ K 2d -« 
+ e^W^Woce K d ' a E d K d ), 



where we use the values of ||//jv||oo,£> R> K > m an d E at the time to- 

Of course if any of the above quantity is infinite then the result is obvious. 
This lemma could appear stupid since we control E(t ) by itself (and with a 
power larger than 1 in addition). But the point is that except for the first 
term, the other two are very small because of the e in front of them so that 
they almost do not count. 

Proof. For any index i, we bound for any ti less than t — e or t 1 = if t < e 

i pti+e 

P = - \E(X t (s))\ds, 
£ Jtl 

which will give the desired result by taking the supremum over all i. So now 
let us fix ti and i (we choose i — 1 for simplicity). 

We introduce the following decomposition among the particles: We define 



Cfc 



3eK(t )2 k ~ L < |Xi(fi) -Xi(fi)| < 3eK(t )2 k \. (2.1) 



Therefore the index k varies from 1 to k = (ln(i?/ '4 e K(to)))/ In 2. And we 
denote by C the rest, that is the set of indices % such that \Xi(ti) —Xi(ti)\ < 
3eK(t ). 

Consequently we decompose the term I 1 into two parts, I± + Ic with 

and Ic is the sum of the same terms over the particles of C . 

Step 1: Stability of the Ck- Given their definition, the Ck enjoy the following 
property, for any i G Ck with k > 1, we have for any t G [t 1 , t ] 



\X 1 (t)-X i (t)\>eK(t )2 



k-l 
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Indeed, we of course know that 



-( Xi (t) - X 1 (t)) 



\V i {t)-V 1 (t)\<2K{t ), 



and then 



\X!(t) - Xi(t)\ > \Xtih) - Xi(h)\ -2 (t - h) K{t ) 
> 3eK(t )2 k - 1 -2eK(t ), 

with the corresponding result since k > 1. Of course the same argument also 
shows that if i G Co then for any t G [ti, to], 

{X^-X^l^beKito). 

Step 2: Control of I±. Using the result from the previous step, we deduce 
that for any % G Cf. with K > 1, 



< 



C2 



-ak 



\Xi(t) - X^t)]" - e a (K(to)y 

On the other hand, we have of course \C^\ < N and moreover \C^\ < 
C s~ 2d K 2d e d 2 dk x ||/ijv||oo,£ according to the very definition of this discrete 
L°° norm (jl.l3j) . Consequently for any a' < d, since e 2d = C/N, interpolat- 
ing between these two values, we get 

\C„\ < CN{K{to)) 2a ' e a '2 a ' k x \\f, N (t , .)||^. 

Using these last two bounds in (|2.2|) and summing up, we obtain 



h<J2 \ Ck \ x N ~* (K(to))~ a e- a 2-° k 



k=l 



ko 



< c mN\\j, E K e 2^ 2 



k=i 



Eventually for any a < a' < d, we deduce that 

h<C \\^ N \\ a J d K 2a '- a e a '- a 2 {a '- a)ko < C \\ii N \\ a J d R a '- a K c 



(2.3) 
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the values being taken at to, which gives the first term in Lemma \2. II Before 
dealing with the remaining term, we point out that here we have never used 
the condition a < 1 and that for this term the same computation would be 
valid for any a < 3. 

Step 3: Redecomposition of Co. For the force induced by the particles in Co, 
we divided again this set into several parts. Set 



Qi = <jeC 



ZeE{t Q )2 1 - 1 < {V^-Vjih)] < 3eE(t )2 l j, (2.4) 

for / > 1. Remark that Q l = if Z > l Q = \n(K(t )/(eE(t )))/ \n2. Therefore 
the rest Qo is defined by 



Qo = <jeC { 



{V^-V^l <3e £(*„)}. 



As before we decompose Iq in a sum of I 2 and a remainder Iq with 

h = Y.T.-J N \ X ^- Xl {t)y (2 - 5) 

1=1 j£Q t J *l 1 3 W IV 1 1 

and for Iq the same sum but on the indices j G Qo of course. 

The idea behind this new decomposition is that although the particles in Qi 

witth I > 1 are close to Xi, their speed is different from V\. So even if they 

come very close to X\ they will stay close only for a very short time. Since 

the singularity of the potential is not too high, we will be able to bound the 

force. 

Step 4-' Stability of the Qi. Just as for the Ck, we may prove that for any 
time t in [t\, t ] and any j G A\ with / > 1 

\V j (t)-V l (t)\ >eE(to)2 l -\ 

This is again due to the fact that 

|V#)-^-(*i)l < f ° \E(Xj(s))\ds<eE(t ), 



so that in fact the result is even more precise in the sense that the relative 
velocity Vj(t) — Vi(t) remains close to Vj(ti) — Vi(ti) up to exactly e E(to). 
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Step 5: Control of I 2 . Given this previous point, for any j G Qi with / > 
and any t G [t%, t 2 ], we have, denoting by t m the time in the interval [ti, t ] 
where \Xj(t) — Xi(t)\ is minimal 



iX^t) - Xjit)] > 
Then, 



\Xi(t m ) — Xj(t m ) \ — ~(t — t m )\Vi(t m ) — Vj(t m )\ 



e 



< C£- 2Q (E(t )) _a 2- a '. 



Summing up on I, we obtain 



/ 2 | <c'£|Q,|ie- 2 «(£;(i„))-«2-"'. 



i=l 

We bound \Q t \ by |Q { | < C WnW^ (K(t Q ) e) d (2 l E(t ) e) d using again the 
definition of the discrete L°° norm. It gives us 

h 

\h\ < C(K(to)) d (^o)) d - a ^ 2a ||^||oo, £ x^2( d - a ) / 

1=2 

d f^U \\d-a\\.. n _2d-2a ( X(tp ) 



< C(K{t,)f(E{t )f- a \\^ N \\^ e e 

< C(K(t )) 2d - a \M\oo,ee d ~ a , 
which is indeed the second term in Lemma \2 .11 



\E(to)e 



Step 6: Control on Iq . The first point to note is that for any j G Qo and 
any t G [t\, to] by the definition (jl.llj) of AE and the stability of Co 

\V 3 (t) - V x {t) - Vjih) - ^(tOI < 5e 2 K(t ) AE(t ). 

Consequently it is logical to decompose (again) Q in Q' U Qq and Iq in 



Qo = U e Qo 



|V#i)- Vi(ti)| >6£ 2 ^(t )AE(t ) 
11 



and Iq> , Iq» the sums on the corresponding indices. 

Then for any j G Q' Q , the same computation as in the fifth step, shows that 



e 



But of course the cardinal of Q' is bounded by the one of Qo and using as 
always the discrete L°° bound 

|g;i<c(K(t )) d (^(to)) d ll/i7vlloo, £ . 

Eventually that gives 

Iq> < Ce 2d ~ 3a (K(t )) d - a (E(t )) d \\fi N (t , .)|U, £ , 

since AE(t ) comes with a negative exponent and being non decreasing, may 
be bounded by the initial value. So Iq' q corresponds to the third term in the 
lemma. 

Let us conclude the proof with the bound on Of course if j G Qq then 
for any t G [t x , t ), 

|^-(*)-ti(t)| < lle 2 K(t )AE(t Q ). 

Now we use the definition (jl.iup of m and the assumption in the lemma to 
deduce that 

\Xi(t) -Xx(f)| > -j- - \V 3 {t) - Vl<*)l > e 2 K(t )AE(t ). 
m(to) 

We bound \Qq \ by Qo which is the best we can do since the discrete L°° norm 
cannot see the scales smaller than e and we obtain 



IqS < Ce 2d - 2a (K(t )) d ~ a (E(t )) d \\v N (to 



i •/ II oo,ei 



which is dominated by the previous term and the third term in the lemma. 
Before ending the proof we wish to note that the condition a < 1 was only 
used to get the second term in the fifth step and in the last step and the 
bound on m was only required in this last step. □ 
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2.2 Estimate on AE 

We may show the following with the same remarks as for Lemma \2. 11 
Lemma 2.2. For any a' with a < a' < 3, assume that 

m{to) < 



12 eK(t )AE(t y 
then the following inequality holds 

AE(t )<C(\\fi N \\£^ d K 1+a ' R*-<* + e d - a -e\\n N \\o 0>£ K 2d - a 
+ e 2d - 3a -P ||mjv|U, £ K d - a E d K d ), 

where we use the values of ||/^jv||oo,e; R, K, m and E at the time t . 

Proof. The proof follows the same procedure as for Lemma \2. II with exactly 
the same decompositions. We have to bound, since as before the choice of 
the indices does not matter 

AJ = 1 lEjX.it)) - E(X 2 (t))\ ^ 



ej tl eP + lX^-XiWl 

We introduce the same decomposition as for the proof of Lemma 12.11 ex- 
cept that now we have two decompositions: One around X\ denoted by a 
superscript 1 and another one around X 2 . So we set for 7 = 1, 2 



3 e K{t Q ) 2 k ~ v < \Xifa) - X^h) I < 3e K(t ) 2 k , 



with Cq and C 2 the corresponding remaining indices. We also denote Co 

cl u c 2 . 

Then of course 

1 r to dt 

ai = £ - to - *> - nx, - x t) \ x ee + w m + £ 

= A/ 1 + A/ Co . 
The second term is easier to bound as we simply write 

Al2 "£^^ (iV \Xi(t) - X 1 (t)\- + N \Xi(t) - X 2 (t)\°) dt 
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Now we do exactly what we did for Ic in the proof of Lemma ETT1 and we get 
as a bound the two last terms in the estimate for E in this lemma, divided 
by s^, which are exactly the two last term in the estimate for AE in Lemma 
12.21 Note by the way that the terms where we sum on Cq for instance but 
with Xi — X 2 in the denominator are in fact even easier to handle. 

Now for Aii, we observe that for i ^ Cq then for any t 
\E(X 1 (t) - Xi{t))-F{X 2 (t) - Xi{t))\ < C|Xi(0 - X 2 {t)\ 

X {n iX^t) - X t (t)\ a+1 + N ^{t) - XW+i) ' 

since it is always possible to find a regular path x t (s) of length less than 
2\X 1 (t)-X 2 (t)\ such that x t (0) = Xi(t), x t {l) = X 2 (t) and \x t (s) -Xi(t) \ is 
always larger than the minimum between \Xi(t) —Xi(t) \ and \X 2 (t) — 
Therefore we are led to make exactly the same computation as for i x with 
a + 1 instead of a which gives the desired result. This is possible only because 
in the estimate on ii, we never use the condition a < 1. □ 



2.3 Control on m and K 

We prove the 

Lemma 2.3. Assume that 
then we also have that 

m(t) < m(0) X e Ct+CeAE(t)+CtiAE(s) ds ^ 

and we may eliminate the e AE(t) term ift > e. 

Note that we still need an assumption on m but it is a bit different (and 
somewhat "harder" to satisfy) than the corresponding one for Lemmas 12.11 
andQ 
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Proof. We consider any two indices i ^ j. Then we write 
d 



e m(s) - Vj( s )\ + lEjXjjs)) - ^(X,( g ))|) 

Since m(t) < £ 1_/3 , the same is true of m(s) and at least one of the quantities 
\Xi(s) — Xj(s)\ and \Vi(s) — Vj(s)\ is larger than therefore 

d ( e \ C e 



ds WXiW-XjW + iViW-vMJ ~ M*)-*i(*)l + M*)-^0OI 

i \E(X<(8)) - E&jsW 



£ /3 + \ Xl (s) - X^s 
But by the definition of AE, see (11.11)1 . we know that for t > e 

J, e» + \X,(s) - X,-(<0| y 
and of course for t < e 



eP + \X i {s)-X j {s) 
Hence, integrating in time, we find 



< 



{XiW-XjW + \v i ( s )-v j ( s )\ - |^(o)— ^-(0)1 +jv;(o)-vs-(o)| 

x e Ct+Ce AE(t)+C f* AE(s) ds 

wich after taking the supremum in % and j is precisely the lemma. □ 

As for K, using the equation that Vi(t) = Ei(Xi(t)), we may prove by the 
same method which we do not repeat, the result 

Lemma 2.4. We have that for any t 

K{t) < K(0) + Ct + CeE{t) + C [ E{s) ds. 

Jo 
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2.4 Conclusion on the proof of Theorem 11.11 

Here (but only in this subsection) for a question of clarity, we keep the 
notation C for the constants appearing in Lemmas 12.11 12.21 12.31 and 12.41 and 

we denote by C any other constant depending only on R(0), K(0) and m(0). 
We assume that on a time interval [0, T], we have (for a given a') 



which we may always do since all these quantities are continuous in time 
(although they may a priori increase very fast). 

Then we show that if T is too small we have in fact the same inequalities 
but with a 3/2 constant instead of 2. By contradiction this of course shows 
that we can bound T from below in terms of only R(0), K(0) and m(0) and 
it proves Theorem II. II with c = C x 2 Sa ~ a+1 . 

First of all, we note that since m(t) < 2m(0), we may apply Lemmas 12.11 
12.21 and 12.31 Furthermore we immediately know from (|1.14j) that 



Let us start with Lemma \2.1\ using the assumption (|2.fij) we deduce that for 
any t G [0, T], 

E(t) < C2 8a '- a (m(0)) 2a ' (K(0)) a ' (R(0)) a '- a + C e d - a + Ce 2d - 3a . 
For e small enough this proves that 



m(t) < 2 m(0), E(t) < 2C 2 8a '~ a (m(0)) 2a ' (K(0)) a> (R(0)) 
K(t) < 2(1 + ^(0)), R(0) < 2(1 + ^(0)), VtG [0, T], 



(2.6) 



HMMIk £ < (8m(0)) 



S(t)< 




8a'- 



a 



(m(0)) 2a ' (ir(0)) Q ' (i2(0)) 



which is the first point. 

Next applying Lemma l2~2l we deduce that for any t G [0, T] 



AE(t) < C. 



From Lemma \2. 31 we obtain that 



m{t) < m(0) x e 1 
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so if T is such that CT < ln(3/2) then we get 

3 

m(t) < — m(0). 

Lemma f2.4l implies that for t G [0, T] 

K(t) < K(0) + CT, 
so that again for T small enough 

K(t)< 3 -(1 + K(0)). 

Eventually thanks to relation ()1.6|) . we know that for t G [0, T] 

R(t) < R(0) + T K{t) < R(0) + CT, 

hence the corresponding estimate for i? provided CT < 3/2. 
In conclusion we have shown that if ()2.6|) holds and if T is smaller than a 
given time depending only on R(0), K(0) and m(0) then the same inequalities 
are true with 3/2 instead of 2. By the continuity of R, K, m and E this has 
for consequence that (|2.fij) is indeed valid at least on this time interval thus 
proving Theorem ll.il 



3 Preservation of 




From the form of the estimate on m in Lemma 12.31 it is clear that with this 
estimate we will never get a result for a long time. Indeed, even assuming 
that we have bounded before K and R, we would have the equivalent of 

rh < m x AE < Cm x m 2+2a ' . 

On the other hand this is somewhat strange since, in the limit, the L°° norm 
is conserved. And this preservation is very useful in the proof of the existence 
and uniqueness of the solution of the Vlasov equation, see for instance [T3] . 
But, how to obtain the analog of this in the discrete case? At this time, we 
just have a bound on ||/z/v||oo,e on a small time, and the bound is too huge to 
allows us to prove convergence results for long time. Of course, this norm is 
not preserved at all because we are looking at the scheme at the scale of the 
the discretization. And in our calculation we do not use the fact that the 
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flow is divergence free, a property that is the key for the preservation of the 
L°° norm. 

So what else can we do? One of the solution is to look at a scale 77 > e, 
with e/rj going to zero as e goes to zero. At this scale, we have many more 
particles in a cell and we will be able to obtain the asymptotical preservation 
of this norm. This will be very useful because it will allow us to sharpen our 
estimate on E and 5E. And with this we will obtain long time convergence 
results. 

Now, we will try to give roughly the idea of the proof in dimension 1 before 
beginning the genuine calculations. We choose a time t and a particle i, and 
look at the square of size 77 in the phase space centered on the particle i. We 
called it S t = {(x,v)\\x — Xi(t)\ < r], \v — V(t)\ < 77}. We want an estimation 
of the number of particles in this square at the time t. For this, we first 
wanted to know where were these particles at the time t — e. During the 
interval of time [t — e, t], the particle j has moved of 

Xj(t) -Xj(t-e) « eVjit) (3.1) 
Vj(t) - Vj(t-e) « eE(t,Xj(t)) (3.2) 

To make it correct, we will have to replace the second left hand side by 
an average on an interval of size e in time, as usual, but we keep this for 
explanation. If we take the values in the right hand side for exact, the 
particles we are interested in were in the following domain S t - e . 

St-e = {{x,v)\\x — Xi(t) — ev\ < r] and \v — Vi(t) — eE(t,x — ev)\ < 77} 

There is two steps to obtain S t - e from St. First, translate for a fixed v, the 
set S t H {v = v} of v, for all v. We call the resulting set S' t . Then, for a fixed 
x, translate fl {x = x} of —E(x), for all x. With these steps, we see that 
S t and S t - £ have the same volume (see figure HJ). But, if we keep the set S' t , 
and do this another time, we will obtain a set with a more strange shape and 
so one. So, we will approximate St- £ by a parallelogram and iterate this step 
from the parallelogram. We will choose the parallelogram 

St-e = {( x , v) \ \x - Xi(t) -ev\ <7] 

and \v -Vi(t) - e(x - XiityVE^Xiit)) - eE(Xi(t))\ < 77}, 

(3.3) 

where E e is an approximation at e of the field, defined by: 

F ( T ) - 1 V x-Xj(t) 
£{X) N^dx-X^l + e)^)- 
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t— ep 



Figure 1: Evolution of S t . 
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We use this approximation to obtain a usable value of V E. This set has 
almost the same volume that S t . And if we begin with a parallelogram, we 
still obtain a parallelogram. What we have to check is that this approxi- 
mation make sense, that means that we forgot or added a non relevant set 
of particles. For this, we need that the parallelogram do not become too 
stretched, in other words that the two sides do not become parallel on the 
figure in dimension one. Remember that we are not exactly interested by 
the volume of the set So (the set that we obtain iterating the process till we 
reached the initial time), but by the number of vertices of the initial net- 
working eZ 2d inside the parallelogram. The volume of the parallelogram is a 
good approximation of this if its width is big with respect to e. 
We avoid this if the slope of the sides of the parallelogram are distinct. 
This is ensured if t and s^2, k<t , £ \'VE £ (kt/e,Xi(kt/s))\ are smaller than 1/2, 
because this two quantities are respectively the tangent and the cotangent 
of the angle between the sides of the parallelogram and the rr-axis. Remark 
that the sum is bounded by tAE(t). 

This limits us in time, so we will obtain the conservation of the ||/ijv||oo,»/ only 
on a short time. But this is not a problem. Call this time T'. After this, 
we choose a new scale rj' sufficiently large, by instance ^ffj. Since we have 
proved the asymptotic preservation of H/Awlloo,??, we can do the same think 
again, replacing e by rj and rj by rj' . And obtain the preservation of ||/ijv|| 00,77' 
on a new small interval of time and so on. 

Now, this is time to do the calculation for the proof of what we claimed 
before. We just need to remark that we draw our pictures in dimension one. 
In this case, we deal with true parallelogram. But we will only prove result 
for d > 2 and in this case we will not use true parallelograms. We will deal 
with sets defined by 



where M is a matrix of A^2d(^), and the norm || • || is defined by 

\\(x,v)\\ = max(|x|, \v\) 

For convenience, we will often decompose the matrix M in four blocks like 
below 





20 



Proposition 3.1. Choose (X , V ) G R 2dn , and A,B,C,D four square matri- 
ces in A4 n {R) satisfying the following conditions that will be called the norm 
conditions in the rest of the paper 

max(||A- Jd||, \\D-Id\\) < 1/2 and max(||C||, \\D\\) < 1/2, (3.4) 

with the dual norm on the matrix. We define M as above and denote by N t 
the number of particles in 



(x,v)\ 



M 



x - X 
v-V 



< rj 



We assume as in all the preceding results, that 

m(to) - mh-mm 

Then, there exist a constant C = C(R,K,E, \\[i\\oo,e) , there exists a position 
X' a speed V , four matrices A' ,B' ,C ,D' and then a matrix M' defined with 
this four blocks such that if we denote iV t _ e the number of particles in 



S; 



t-E 



[x,v)\ 



M' ■ 



x-X' Q 



< V + Ce{r ] fS + e) 



we have the following inequality 

N t < N t _ £ . 

Moreover, we have 

max(||£> - D'||, \\D - D'\\, \\C - C% \\D - D'\\) < Ce. 
Proof. We divide it in two steps. 

Step 1: Estimate on Vj(t — e) — v — eE(Xj(t)). Choose j G {1, . . . , iV} such 
that Xj(t) G S t . We will first work on the velocities and then integrate our 
estimation to obtain what we need on the positions. We have 

Vj(t-e) -Vj(t) = el E{Xj{t-se))ds 

Jo 

= el E(Xj(t- se)) - E £ (Xj(t- se))ds 
Jo 

+e / E e (Xj(t- se))ds 
Jo 

= 1 + 11. 
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We need to bound the first term I. The approximation error is 



EiXAt - se)) - EJXAt - se)) = — Vfi — -- — 1 , XM — 
v A ;; eV n ;; N^KlXjit) -X k {t')\ 1+a 



' \XAt')-X k {t')). 



(iX^-X^l + e 



l+a I v J v 



To compute this term, we use again the same decomposition of the phase 
space as in lemma l2~2l If k is such that \X k — Xj\ > e, we can bound a 
difference in the sum by 2 1+a e/ {\X^ — Xj\) 1+a and then compute it dividing 
the phase space in diadic subset. We obtain that 

2 1+a e/(\X k -X,\r<Ce. 

\X k -X \>e 

The constant C that we obtain here depends of R, K, E, H^tvII oo,e exactly as 
in lemma 12721 For the others terms, those were \X k — Xj\ < e, we bound the 
difference by the sum of the norm of the two term and do the same estimation 
as in ()2.2|) . It is possible because we do a mean over a interval of time of size 
e. We obtain 

|/| < Ce 2 + C'e d+1 - a . 

again with the same dependence for C . In the following, we will assume 
to simplify the presentation that d > 2. In that case, we may replace it by 
/ < Ce 2 and this bound will be enough. But all is still true if d = 1. In that 
case, we replace the last estimate by I < Ce 2 ~ alpha . Since — a > 1, all our 
estimates will remain valid. 

Now, we approximate the second term 77 by the same expression where we 
replace E e by its first order linearization near Xq. We defined E l £ m by 

E l J n (t,x) = E e (t,X ) - VE £ (t,X ) ■ (x-X ). 

We have to estimate the difference between E £ and E l e m . For this, we need 
to bound a sum of terms of the following type 

X — Xk Xq — Xfc 



(\x-X k \ + ey+« (\x -x k \+ey+« 

x - X k 



(Ix-Xkl + s) 1 - 



{t,X )-{x-X ). (3.5) 
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For this, we find a path J(s) between x and X (in others words 1(0) = x 
and 1(1) = Xq) of length smaller than 4\x — X \ and such that \I(t) — X k \ > 
min(|x — Xk\, \Xq — X k \). The term to estimate may be rewritten 



1 v Qi^V )^)) 

x - X k 



(t,X )) ■I'(s)ds, 



X\x~x k \ + ey+ a , 

and each difference (without the multiplication by I'(s)) may be bounded by 

C 

where C is just a numerical constant and also by 

C 





I(t)- 


-x \ 




min( \x — 




1 


x - 


-x k \)^ 





x - 


- x k \ 


1+7 


min( \x — X k \ 




x 


-Xkiy+a+'r 



mm(\x-X k \,\X -X k \y+«' 

where the constant C is again independent of the problem. Doing a classical 
interpolation between this two terms, we may bound it for every 7 6 (0, 1) 
by 

\I(t)-X Q \^ 
min(\x-X k \,\X -X k \y+^- 
Thus, we may bound the approximation error by 

-E 

N ^ 

Now, we do the same computation as in the estimation of AE. We obtain, as 
in lemma (J2.2|) . that for every 7 = (3 — 1, there exists a constant C depending 
always on R, K, E, \\/j,n\\oo,s such that 

I- / E £ (s, Xj(t - se)) - E l E in (s, X 3 (t - se)) ds\ < C(\X 3 j(t) - X f + e). 

£ Jt-e 

We approximate the term // by 

\II-e [ Ef n (s,X,j(t - se))\ds < Ce(\Xj(t - se) - Xof + e). 
Jo 
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Moreover, we know that the approximated force field is Lipschitz with a 
constant depending on R, K, E and ||//jv||oo,e 011 our small interval of time. 
This can be seen just by doing the same calculation that gives the bound on 
AE. Indeed, it is the same because regularization at order e or estimation 
on AE at order e are similar. Moreover, we may replace the value of E l ™ 
taken at Xj(t — se) by those taken at Xj(t). Because of the bound on the 
speed, this will only introduce an error of order e. Actually, if we denote 
E = l/e ft_ £ E £ (s,X )ds and VE = 1/e f^_ £ VE £ (s, X ) ds, we have the 
following inequalities 

\V j (t-e)-V j (t)-e(E + VE-X j (t-e))\ < Ce{\Xj{t) - X \ + e), (3.6) 
\Xj(t - e) - Xj(t) - eVj(t - e)\ < Ce 2 , (3.7) 

where the constant C has the same dependance as before. 

Step 2: The new parallelogram. Now, we obtain an approximated parallelo- 
gram that contains at time t — e almost all the particles that are in St at time 
t. We look first at the two first blocks of the matrix and define the following 
linear mappings 

A' = A + eB ■ V£ e , 
B' = B + sA. 
The center of our new parallelogram is given by 

X' = X - eV , Vq =V - f E £ {s, X ) ds. 

Jt-e 

The previous calculation tells us that 

\A' ■ (Xj(t - e) - X' Q ) - B' ■ (Vj(t - e) - Vfi\ < rj + Ce(\Xj(t -e)-X^ + e). 

Of course we obtain the same for the second line, and if we define also 
C = C + eD ■ VE £ , D' = D + eC and the associated M', we obtain 

< Ce(\\Xj(t) - X ,Vj(t) -V f + e). 

To conclude, we just have to prove that \Xj(t) — X \ is of order rj. But, this 
is true because our initial parallelogram is not too stretched thanks to the 
conditions on the norm of A, B, C, D. The following lemma is more precise. 



M >( X j {t-e)-X' \ 
I V 3 (t-e)-V> ) 
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Lemma 3.1. Assume that the four matrices A, B, C, D satisfy max(||A — 
Id\\, \\D — Id\\, \\C\\, \\D\\) < 1/2. We use here the dual norm on matrices, 
\\A\\ = sup|x|<i \ A ■ X\/\X\. We denote as usual by M the matrix composed 
of the four blocks A, B, C, D, and 

S = {(x,v)\\\M ■ {x,vf\\ < 77}. 

Then 

S C 5(0,2??). 

Proof of the lemma. Choose an x G S. With the definition of our norm || ■ || 
it means that 

\Ax + Bv\ < 77, 
\D v + C x\ < rj. 



This implies 



\Ax\ — \Bv\ < 77, 
\D v\ — \C x\ < rj. 



With the assumption on the matrix, we obtained 

N-^(N + H) < V, 
\v\ ~ + \v\) < 77. 



We add this two lines and obtain 

\x\ + \v\ < 2i]. 

□ 

There only remains to prove the estimate on the determinant of M' . But, 



M'— 



A + eBVE B + eA 
C + eD-VE D + eC 
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Then, 



A+(*/fi\ aJ A + sBVE B + eA\ 

d6t(M) = det {c + eDVE D + eC) (3 ' 8) 

, / A-e 2 AVE B + eA \ 
- aet {c-e 2 CVE D + eC ) 

(3.10) 

To obtain the second line, we substract the second column multiplied by VE 
to the first. Then, we see that the new determinant is e 2 close from the 
preceding. This gives us the expected bound on the determinant of M' . □ 

With the help of the proposition, we prove the almost preservation of the 
L°° norm at the scale rj, which is stated in the theorem below. 

Theorem 3.1. There exist a time T' , a constant C both depending on R(T'), 
K(T'), E(T'), supt<T'\\fJ>N{t)\\oo t e) suc h that if t < X" and 

m(T') < ' 



12 e K(T') AE(T') 
the following inequality holds: 

Wf^N ^ J 1 1 oo ,77 < Mint + C{rf + e/rj). 



Proof. Let us choose a box S t = {(x, v) \ \\x — X ,v — Vq|| < rj}. We can 
find a parallelogram containing at time t — e the particles that are in St at 
time t. And we can iterate this process till the parallelogram do not become 
too much stretched. In others words, till the conditions 13.41 are satisfied by 
the matrix M. As at each step, the matrix move from at worst C'e, we may 
iterate the process on an interval of time of length 1/2C". So we obtain a 
parallelogram So not too stretched, if t < 1/2C. We made ([t/e] + 1) steps 
to reach the time 0. The last is necessary of length less than e, but as in the 
proof of the theorem one, this raises no difficulty So, we know that Sq is not 
too stretched, but we have to check that this parallelogram So is not too big. 



Step 1: The size of So We define real functions g £ by g £ (r)) = 77 + Ce(e + rf\ 
We introduce it because if S t ' is defined by 



S f , 



[x,v)\ 



M ■ 



x ■ 

v 



Xo 
Vo 



< 
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Then S t >- £ is denned by a similar condition with new M, Xq,Vq where the 
right hand side of the inequality is replaced by g e {rj). And, so So is defined 
by a similar conditions, where the right hand side n is replaced by 

rf = 9e o • • -og^T]), 

where the dots means \t/e] + 1 times. We need to control this quantity We 
define r n = g™{rj) where the exponent means composed n times. This r n 
is the quantity in the right hand side of the definition of St-ne- From the 
formula r n+1 = g(r n ) we expect that r n « r\ + Cne{e + rf). To prove this 
rigourously, we introduced a n defined by r n = n + Cne^e + n 13 ) +a n . Provided 
a n < rj, it satisfies the following relation: 

a„+i < (1 + C'erf- X )a n + C C'e 2 n^-\e + rf), 

where C is a numerical constant. We can bound ne by t since we will only 
iterate the process till n = [t/e] + 1. The previous inequality becomes 

cwi < (l + C'e^- l )a n + CC'e^-\e + ^). 

Since «o = 0, we obtain that 

a n <((l + C'e V p - 1 ) n -l)C(e + ^). 



This gives us 



r{ < T 1 + {l + e c '^~ 1 -l)C{e + r} P ) (3.11) 
< n + C{e + rf), (3.12) 



with a new constant for the last line. 

Step 2: Covering of So by balls of radius e/2. 

Lemma 3.2. let S be a parallelogram defined as above by 



S=\(x,v)\ 



M 



x - X 

V-Vq 



< rj 



with M composed of the four blocks A, B, C , D satisfying the assumption 
\3.4\ and also det(M) < 1 + Ce. Then, there exists a constant C depending 
only on C such that S can be covered by [e~ 2d (Vol(S) + Cerf d ~ 1 )} balls of 
size e. In others words, there exists a finite set P of cardinal less than 
e- 2d (Vol(S) + C'erf d - 1 ) such that S C U peP B(p, e) 
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Proof of the lemma. We define 



S. 



2s 



s 



Is 



(x,v)\ 
(x,v)\ 



M 



M ■ 



x — Xq 

v-V 

X- Xq 

v-V 



<n + 2e 
<n + 4e 



and P = eL fl S^- We look at the set P £ consisting of the union of all the 
balls (for the norm max(|x|, \v |)) of size e centered at points of P, that is 
P £ = P + B(0, e). We will show that this set is included in Sf £ . For this, we 
choose (x, v) G P £ . We associate to this point the couple (m,n) such that 
II (x — em, v — en) II < e. Then, 



M ■ 


C) 


< 


M 











M 



em 
en 



x — em 
v — en 

< \\M\\e + n + 2e 

< n + Ae. 

In the last line, we use ||M|| < 2. This inequalities comes from the condition 
13.41 Therefore we have the inclusion P £ C S£ £ . 

Moreover, if we choose a point (x, v ) G S, we can find a point (em, en)of el? d 
such that || (x — em, v — en) \\ < e. As above, we have 

< n + 2e. 

Thus, e(m, n) G P. That proves that S C P £ . So, we have the inclusions 

S G P £ G S^ £ . 

The first is the recovering we want. The second gives us an estimate on the 
cardinal of P. Comparing the volume of P £ and S£ £ we obtain 



M ■ 


( em \ 




\en ) 



(2e) M \P\ < det(M)(n + 4e) 2d . 
If M satisfies det(M) < 1 + Ce, that gives us 

\P\ <e- 2d n 2d -\n + C'e). 



□ 
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Step 3: Conclusion of the proof. We choose t as in the previous step and a 
box S t in the phase space of size r\. We define the parallelogram So as above. 
It exits and is not too stretched because we choose t <T'. From Step 1, we 
know that 

Vol(So) < (l + Ce)(r ] + C(E + r ] ( ')) 2d (3.13) 
< rf d (l + C{r l l3 + e)). (3.14) 

And since it is not too stretched, we may used lemma I3~2l and cover it by less 
than e~ 2d (Vol(S ) + C'erf d ~ l ) balls of radius e. Actually, that means that 

N t < N < 11^(0)11^^(1 + CV + e), 

and dividing by rf d we obtain 

||A*iv(*)IU,oo < |M0)|| £)O c(l+C(^ + £)). 

□ 



3.1 New estimates on E and AE 

The almost preservation of the ||/xjv||oo,r; norms will enable us to prove a new 
estimate on E. We can obtain it by doing the same separation of the position 
space in dyadic cells, but we begin with cells Cp. satisfying 



Cu 



3 V K(t )2 k ~ L < \Xiih) - X 1 (t 1 )\ <3 V K(t )2 h 



That gives a first term in \\/an 

W^n^R - . Next we decompose Co in a 
union of Ck (the "true" cells with size e); Since the index k goes only up 
to log (77/e) /log 2, the corresponding term is less than ||/xjv||oo,e K 2a ~ a r] a ~ a . 
The remainder term coming for Co is dealt just as in Lemma l2~T1 and it yields 
the two same terms. 
That gives us the following lemma 

Lemma 3.3. For any a' with a < a' < 3, assume that 

m{to) < 



12 eK(t )AE(t o y 
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then the following inequality holds 

E(t ) < C (\\fi N \\ a J d K a ' R a '~ a + \\ f iN\\^K 2a '- a r } a '- a 

i ^d—a II,, || T^2d— a . Jld—Za II,, || isd— a~Ti a 



where we use the values of ||//jv||oo,e> R, K, m and E at the time t . 

The only non-negligable term in this estimate is sub-linear if a 1 is chosen 
sufficiently close to a. 

Of course we can perform the same changes in the proof of AE to get 
Lemma 3.4. For any a' with a < a' < 3, assume that 

1 



12eK(t )AE(t ) 
then the following inequality holds 

AE(t ) < C {\\fi N \\£+ v Q ' )/d K 1+a ' R a '- a + \\v N \\£+ a 'V d K 1+2a '- a r) a '- a 
+ e d ~ a -P |K|U £ K 2d ~ a + e 2d - 3a -P H^IU,, K d ~ a E d ), 

where we use the values of ||//jv||oo,e? R, K, m and E at the time t . 
3.2 Proof of Theorem fT2l 

Let us fix any time T > 0. The aim is to show that we have bounds for R, 

K, E and m, uniform in N on [0, T]. 

Next we choose r/o — s 1 ^ 2 for instance and i]' = e 1 ^ 4 . 

Since for any N the quantities R, K, E and m are continuous in time, we 
may define Tjy < T as the first time t (if it exists) such that one of the 
following inequality at least is not true for some integer M to be chosen after 



r = T(R(t),K(t),E(t),mp\\ f i N \\ 00te ) > 

s<t M 

m ^ ^ io e^aeV^ ' C(R(t),K(t),E(t),suv \M\oo,e) < ^ 1/8A/ , 
12 e K(t) AE(t) s <t 

E d - a (m(t)Y 2d (K(t)) 2d - a < e , e 2d - 3a (m(t))- 2d (E(t)) d (K(t)) d - a < e p . 

(3.15) 
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The quantity T' and C are the time and constant defined in Theorem 13.11 
Therefore on [0, Tjy] all inequalities (|3.15|) are true and we may apply both 
Theorem 3.1 and Lemma f3. 31 

We define tj = %T' and r\i = t/q x r l with r = e -1 / 4 * 7 so that t)m = v'- We 
are going to apply M times Theorem 13.11 once on every interval t;] (of 
size less than T") and with rj = rji and e replaced by rj^i. That gives 

sup \\tiN(t)\\oo,TH < IMk^ +C(E(t i ),AE(t i )) {T]J + e 1/iM ), 
te[ti-i,U] 

and consequently thanks to ()3.15|) 

sup IKWIk,/ < \\I2 N \\^ £ + C(E(T N ) } AE(T N )) M e 1 ^ 1 <2\\fi° h 



t<T N 

(3.16) 

independently of iV (and T^). Now we apply Lemma 13.31 at time T/v and 
because of (13. 15 J) , we obtain 



(3.17) 



N 



E{T N ) < C \\li N {T N )\\t{t {K{T N )Y (R(T N )r'-° 
<C(K(T N )f (i?(T^)) Q '- Q , 

using (|3.16jl . As T N > e, Lemma E31 implies that 

fT N _ __ 

K{T N ) < K{0) + C E(t) dt < K{0) + CT N E{T, 
Jo 

From this inequality, we immediately deduce that 

R(T N ) < R(0) + T N K(0) +CT 2 N E(T N ) <CT + CT 2 E(T, 

Inserting these last two inequalities in ()3.17|) . we find 

E{T N ) <CT + CT 2 (E{T N )) 2a '- a . 

Since 2a' — a < 1, there exists a constant C{T) depending only on T and 
the initial distribution such that 

E(T N ) < C(T), K(T N ) < C(T), R(T N ) < C(T). (3.18) 

We are almost ready to conclude, we only need to apply once Lemma 13.41 
and by dSHSJ) , and (j3~TH|l 

AE(T N )<C(T). (3.19) 
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Inserting (J3.19)) in Lemma [2.31 we eventually get 



m(T N )<C(T). 



(3.20) 



Together (j3~T£|) . (jHUHD and (j3~27Ijl imply that all the inequalities of (j3~T5jl 
are true with a factor 1/2 at time T/v, provided N and M are large enough. 
Therefore ()3.15|) is still true on at least a short time interval after T/v and 
that means that necessarily T/v = T. The consequence is that (j3.18J) . ()3.19|) 
and ()3.2U|) are true on any time interval [0, T] which is exactly Theorem 11.21 
Finally note that we have implicitly used the short time result when we said 
that T N > e. 

4 Convergence of the density in the approx- 
imation 

The existence of the bound on R, K, E, AE and ||//jv||oo,jj implies the weak 
convergence of the distribution fi N to a weak solution of the Vlasov equation 
and Theorem 11.31 is only a consequence of Theorem 11.21 and the following 
proposition 

Proposition 4.1. Let /i/v be the distributions associated with the solutions to 
(jl.lj) . We assume that the initial conditions fi° N converges weakly in M l (M. 2d ) 
to some fo G L 1 PI L°°(M. 2d ). We choose a time T > 0. Assume furthermore 
that there exists a constant C (T) independent of N such that 



where r\ depends on e and N and goes to zero when e goes to zero. Then, 
yUjv(t) converges weakly to f(t), a solution to the Vlasov equation with initial 
conditions f°. 

Proof. We recall that the distribution of the particles satisfies the Vlasov 
equation in the sense of distribution. Moreover, the sequence /i^ is bounded 
in C([0, T], M 1 (R 3d )). Up to an extraction, we may assume that [in converges 
weakly to some / E L°°([0, T), M 1 (R 2d )). Moreover, the fact that Halloo,,, is 
bounded implies that / G L°°. To see this, we choose a regular test function 
cf) with compact support. We have 



sup( J R(T), K(T), E(T), AE(T), W^J) < +oo , 




i=l 
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Now, we define p v (x,v) = Xci^/v^/v) where xc is the characteristic func- 
tion of the set C = {(x,v)\\\(x,v)\\ < 1} and we write 

1 N 

1=1 

1 N 



+ n w)) - $ * ^ (*)))• 



1=1 



The first term is J <p*p v (x, v) dpN{x, v ) = / 4>{pN*Pri) dxdv. So it is bounded 
by ||0||i||0A t Ar*Pr?||oo- But \\<Ppn* P v \\oo is exactly ||/ijv||oo,7/- The second term 
is easily bounded by ?y||V$||oo- Putting all together, we obtain that 

(A*JV,$) < ||Piv||oo,r,||$||l + ^||V$||oo- 

At the limit, 

(/,$) < liminf H/XArlloo^H^lli, 

which means that / e and that ||/||oo < liminfTv^oo ||a*jv||oo,tj- 

The passage to the limit in the linear part of the equation does not raise any 
difficulty For the term in F • V„/, we need a strong convergence in the force. 
We denote by the force induced by / and by F N the force induced by p^ 



/x — y 
\x - y\ L+a 
1 N 



x - Xj(t) 
\ x - Xt (t)\ 1+a ' 



We have 



i r< 



F w (X,(t))-F 00 (X l (t)) dt = h + I 2 + h 

y-x t (s) 



i r +£ /• 

e Jt J\y- 



X i{ s)\>r \V ~ X^S)^ 1 



d{pN - f)(y) ds 



V ' XM df ( y)ds, 



i r t+£ r 

£ Jt J\y 



. Xi(mr \y - XiW+i 
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for all r > 0. The first term 1% in the right hand side always goes to 
zero because fi^ converges weakly to /. The second term is dominated by 
ll/llco Jbior) dy/\y\ a , a quantity which is less than CH/HooT^ - ". The last one 
is the field created by the close particles in the discrete case. To estimate it 
we will use the lemma l2~H replacing R(T) by r. This gives 

+ e 2d - 3a \\ f i N \\ OD , £ K d - a E d K d ) < Cr a '~ a . 

And these bounds are independent of N or i. 
Then, letting e going to and then r, we find that 

1 f t+£ 

sup-/ \F N (Xi(s)) - FcciXiis))] ds^O as e -> 0. (4.1) 

i,t £ Jt 

With this strong convergence, we are able to prove the convergence of the 
term F N ■ V v fiN ds towards F^ ■ V v f in the sense of distributions. We choose 
a test smooth test function <p with compact support and compute 

J— ( / F oa (t,x)-V v (f>(t,x,v)f(t,x,v)dxdv 

JO ^ Jx,v 

N 



- Ys F N(t,X t (t),V t (t)) ■ V v 0(t, V^))) dt (4.2) 



We separate J in J\ + J 2 , with 



J\ = / / -Foofoz) -V v (j)(t,x,v)d(f - fi N )(,x,v)dt, 

JO Jx,v 



and 

„t N 



h = f ^F^t^X^^it^-FN^X^Mt^-^M^X^^^dt 
»=l 

Because of the continuity of i 7 ^, vanishes as e goes to zero. To show 
that J2 vanishes as well, we decompose it in M = [T/e] + 1 integrals on M 
intervals of time with length e. The last interval is of length less than e, but 
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that does not any difficulty and we do as if it were of length e. We obtain, 

M p(k+l)e N 

k=l Jke i=l 

•V„0(*, ^(t)))dt (4.3) 

M Mk+l)e N 
k=l ^ ke i=l 

This sum may be bounded by 

CT sup - / + \F N (Xi(s)) -F^X^ds, 
i,t £ Jt 

a quantity which goes to zero according to (|4.1|) . Thus, J goes to zero when 
£ goes to zero and the proof is done. □ 
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